Momentum dependence of orbital excitations in Mott-insulating titanates 
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High-resolution resonant inelastic x-ray scattering has been used to determine the momentum de- 
pendence of orbital excitations in Mott-insulating LaTiOs and YTiOs over a wide range of the 
Brillouin zone. The data are compared to calculations in the framework of lattice-driven and 
superexchange-driven orbital ordering models. A superexchange model in which the experimentally 
observed modes are attributed to two-orbiton excitations yields the best description of the data. 
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Valence electrons in transition metal oxides exhibit 
a large variety of ordering phenomena that are associ- 
ated with unusual macroscopic properties [![. In or- 
der to elucidate the origin of these properties, experi- 
mental research beginning in the 1950's has systemati- 
cally investigated the relationship between the occupa- 
tion of d-orbitals on metal ions and the magnetic or- 
dering pattern. The leading paradigm for the inter- 
pretation of these data is the "Goodenough-Kanamori" 
framework, which is based on the assumptions that the 
orbital occupation is determined by electron-lattice in- 
teractions alone, and that the relative orientation of or- 
bitals on neighboring sites controls the exchange inter- 
actions between magnetic ions While experiments 
on most transition metal oxides are in accordance with 
this paradigm, it has been called into question by recent 
data on the Mott insulators LaTiOa and YTiOa with the 
seemingly simple electron configuration 3d}. Not only 
does the spin dynamics of these compounds appear to 
defy a quantitative description in terms of the standard 
lattice-driven orbital ordering scenario 0, 0] , thermody- 
namic data Q on LaTiOa have even been interpreted 
as evidence of a novel spin-orbital liquid state [5| that 
is qualitatively incompatible with this scheme. An al- 
ternative approach that emphasizes the many-body su- 
perexchange coupling between spin and orbital degrees of 
freedom while treating the local electron-lattice interac- 
tion as a perturbation, can explain some of these experi- 
mental findings 0,0]. Other experiments, however, have 
been interpreted in terms of the rigid orbital ordering 
pattern predicted by the Goodenough-Kanamori frame- 
work 0,11, 0, [Uj m, ■ Further experimental work 
is therefore required in order to discriminate between the 
competing theoretical scenarios and to develop a compre- 
hensive understanding of the electronic structure of these 
prototypical orbitally degenerate Mott insulators. 



Momentum-resolved spectroscopies, which associate 
specific energy scales directly with the chemical bond- 
ing pattern, are among the most powerful experimental 
probes of transition metal oxides. In particular, the mo- 
mentum q dependence of low-energy (~ 10 meV) spin 
excitations determined by neutron spectroscopy 0, 0| 
and high-energy (~ 1 eV) charge excitations determined 
by electron 14 1 , photoemission , and x-ray spec- 
troscopies has yielded deep insights into the electronic 
structure of Mott insulators. Recently, orbital excita- 
tions with intermediate energies (~ 250 meV) have been 
detected by Raman 17 1 and infrared [ll| spectroscopies 
in Mott-insulating titanates, but these methods are lim- 
ited to g = 0. Because of the much larger photon wave 
vectors involved in the scattering process, resonant in- 
elastic x-ray scattering (RIXS) can in principle serve as 
a probe of the momentum dependence of the orbital ex- 
citations [ii, 0, HH. This is particularly informative 
for the titanates, because of sharply divergent predic- 
tions by the two theoretical scenarios. Whereas orbital 
excitations in the lattice-driven theory are weakly dis- 
persive because their energies are dominated by the lo- 
cal Jahn- Teller coupling [3, [l^, the superexchange 
model predicts collective modes with a substantial dis- 
persion Instrumental constraints have, however, thus 
far precluded RIXS experiments on orbital excitations 
with q 0. 

We have overcome these limitations at the newly de- 
veloped ADRESS beamline at the Swiss Light Source 



23l | and report RIXS data with photon energies near 
the Ti ^2,3 absorption edge of LaTiOa and YTiOs over 
a large fraction of the Brillouin zone. The experiments 
were performed on high-quality YTiOa and LaTiOs sin- 
gle crystals (grown as described in Ref. iS) with magnetic 
transition temperatures Tc = 27 K and Tjv = 146 K, re- 
spectively. The YTiOs (LaTiOs) crystal was untwinned 
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(partly twinned in the ab plane of the orthorhombic space 
group Pbnm) . Since the orthorhombic distortion is small, 
we use the pseudocubic notation, i.e. a' — a/\/2 , b' — 
fo/v^, and c' — c/2. The crystals were cut and polished 
with surface normals along the [110] and [100] directions, 
respectively, for LaTiOs and YTiOa. This allowed inves- 
tigations of orbital excitations with momenta along these 
two directions in the Brillouin zone. The RIXS data were 
taken at the ADRESS beamline, where an energy reso- 
lution of 55 ± 2 meV was obtained using a combination 
of focusing elements and high-precision optical gratings 
(i^ . [2^. X-ray absorption spectra taken to select the 
energy for the RIXS experiment were identical to those 
presented in Ref. (2^. To study the momentum depen- 
dence of the orbital excitations, the incident photon en- 
ergy was adjusted to either 452.0 eV or 454.6 eV at the 
Ti L3 edge (corresponding to excitations from the 2^3/2 
core level into the 3d t2g and eg levels, respectively), and 
the scattering angle 29 was changed in a range between 
50° and 130°. Data were taken at T = 13 K (i.e. in 
the magnetically ordered states) and at room tempera- 
ture. The incident x-ray beam was polarized either hor- 
izontally or vertically (parallel or perpendicular to the 
scattering plane, respectively). While excitations of com- 
parable intensity were detected in both polarization ge- 
ometries, elastic scattering arising from static disorder, 
surface effects, and partly from thermal diffuse scattering 
is strongly enhanced for vertical polarization. The elastic 
peak in vertical polarization was therefore used for cal- 
ibrations of the photon energy (performed periodically 
between inelastic scans) and a determination of the en- 
ergy resolution. The spectra shown below were taken in 
horizontal polarization. They are sums of partial spectra, 
each accumulated for 10 minutes for a total of 4 hours. 

Figure 1 shows the t2g excitation spectra (incident pho- 
ton energy 452.0 eV) of LaTiOa and YTiOa at room tem- 
perature and at T = 13 K in horizontal polarization with 
29 = 90° (panels a,b) . In this scattering geometry the dif- 
fuse elastic peak is strongly suppressed, and the signal is 
dominated by inelastic scattering. All inelastic peaks ap- 
pear on the energy-loss side, as expected because the ex- 
citation energy of ^ 250 meV greatly exceeds ksT. The 
peaks disappear when the energy of the incident beam 
is tuned away from the resonance condition, demonstrat- 
ing that they arise from excitations of the 3(i-electrons of 
the Ti ions. The inelastic peaks sharpen with decreasing 
temperature, and their intensities increase. The same 
effect has also been observed by Raman light scatter- 
ing [13] • Interestingly, the "elastic" peaks at T = 13 K 
are not centered at zero energy, but at 40 (16) meV for 
LaTiOs (YTiOs). The shift closely corresponds to the 
energy range of (anti-)ferromagnetic magnons previously 
detected by neutron scattering [2, [1] . At room tempera- 
ture, where both compounds are paramagnetic, the peak 
centers shift to 22 (0) meV, supporting an interpretation 
in terms of magnons. With further improvement of the 
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FIG. 1: (Color online) RIXS spectra of (a) LaTiOa and (b) 
YTiOs at room temperature and at T = 13 K. The spec- 
tra were taken at an incident energy Ei = 452.0 eV (i.e. at 
the t2g level of the L3 excitation of the Ti-ions) in horizon- 
tal polarization of the incident light at a scattering angle 20 
— 90°. The lines are the results of fits to the Voigt line- 
shape. (c,d) Evolution of the RIXS spectra as a function of 
29. Note that the excitation momentum q = 2ki sin 6, where 
ki is the incident photon wavevector. The data along [100] 
([110]) thus span a range of 11.5 to 25% (8.5 to 17.5%) of 
the diameter of the Brillouin zone. The weak features in the 
spectra arise from statistical fluctuations of the count rate, 
which were smoothed out by the detector read-out algorithm. 



energy resolution, RIXS may thus become widely appli- 
cable as a probe of dispersive magnon excitations. 

The momentum dependence of the RIXS spectra at 
T = 13 K is shown in Fig. lc,d. Kinematic and instru- 
mental constraints restrict the q-range to 18.2% (25.8%) 
of the Brillouin zone for YTiOg along [100] (LaTiOg along 
[110]). Note that the Brillouin zone boundary is at 50%. 
In a first step of a quantitative analysis, the data were fit- 
ted to Voigt profiles, which correspond to a convolution of 
a Lorentzian with the Gaussian spectrometer resolution 
(lines in Fig. la,b). The resulting peak positions and am- 
plitudes of both t2g and eg excitations are plotted in Fig. 



2, along with previously reported [17|, |20| data for q = 0. 
Only a small variation of the peak positions is observed 
in the g-range covered by our experiment (Fig. 2a). This 
behavior (which is very similar in both compounds) is 
incompatible with the strongly dispersive single-orbiton 
excitations predicted by the superexchange model. Our 
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FIG. 2: (Color online) Momentum (q) dependence of (a) the 
energy and (b) the amplitude of the RIXS signal of LaTiOs 
and YTiOs at T = 13 K. g is given as a fraction of the size of 
the Brillouin zone (BZ) along the high-symmetry directions 
quoted in the text. Note that the BZ boundary is at 50%. 
The g = data obtained from Raman light scattering [iJl are 
also shown in (a). The lines are guides-to-the-eye. 

calculations below demonstrate, however, that the RIXS 
intensity of two-orbiton excitations can exceed that of 
single-orbiton modes. The intensity maximum of the 
two-orbiton continuum is determined by a convolution 
of two single-orbiton states with momenta adding to q, 
and its dispersion is expected to average out. 

In order to obtain a quantitative description of the 
contribution of one- and two-orbiton excitations, we now 
consider the scattering processes in detail. We first fo- 
cus on a superexchange model for YTiOs (Ref. @), and 
then compare the results to those of a lattice-driven sce- 
nario. In RIXS, a 2p electron of the Ti-ion is lifted to 
the 3d level, modifying the electronic environment. In 
particular, the core-hole potential effectively reduces the 
Coulomb repulsion U by the amount Uc, and therefore 
modulates the superexchange interaction Jse = 4t^/[/. 
This exchange-bond modulation mechanism was success- 
fully used to explain the observation [24] of two-magnon 
scattering by RIXS [25, 26, 27]. In the titanates, modi- 
fication of Jse causes either one or two orbital flips on 
neighboring Ti~sites. Along the lines of Ref. 121|, the fol- 
lowing effective scattering operator can be derived: 

d, = ^e-^-(jini;) -f J.i^,) , (1) 

i,S 

where Ji = ^ - and J2 = ^ + ^ - ^. The 

operators n and ^ -^+5 were taken from Ref. and con- 
tain the density and pseudospin one-half orbital opera- 
tors whose explicit form depends on the spatial direction 
7 = a, 6, c of the exchange bonds. The first term Ji repre- 
sents single-site processes, whereas the second term acts 
on two sites. In terms of orbitons, the latter term only 
creates two orbitons, while the former can create both 




-0.6 -0.5 -0.4 -0.3 -0,2 -0.1 0.0 0.1-0.6-0,5 -0,4 -0.3 -0.2 -0.1 0.0 0.1 -0.6-0,5 -04 -0.3 -0.2 -0.1 0.0 0.1 



Energy (eV) Energy (eV) Energy (eV) 

FIG. 3: (Color online) Comparison of the model calculations 
with the experimental RIXS data of YTiOa for different scat- 
tering vectors q along the [100]-direction: (a) superexchange 
bond modulation, (b) local shakeup of collective orbital exci- 
tations, and (c) local crystal field excitation model. 

one and two orbitons. The single-orbiton contribution 
scales with the orbital order parameter which is strongly 
reduced @ due to the orbital frustration inherent in the 
superexchange model. Therefore, the two-orbiton re- 
sponse dominates the scattering intensity. Moreover, for 
the [100]-direction the Ji-scattering channel is forbidden 
for symmetry reasons, because the number of j/z-orbitals 
is conserved on each [100]-plane independently. 

The results computed for the bond-modulation RIXS 
process are shown in Fig. 3a. The scattering intensity 
was calculated using the orbital dispersion relation from 
Ref. |6i, with the orbital exchange strength riJsE — 75 
meV (close to the value of 60 meV estimated in Ref. @) . 
In order to facilitate comparison with the experimen- 
tal data, we introduce an orbital damping 7 = 30 meV 
due to orbiton-orbiton (superexchange model) or orbital- 
phonon interactions (Jahn- Teller model). The spectra 
are further convoluted with the resolution function of the 
spectrometer. The lineshapes of the calculated and mea- 
sured RIXS spectra agree very well. 

Figure 3a shows that the bond modulation mechanism 
alone fails to describe the evolution of the absolute inten- 
sity of the RIXS signal with increasing q. However, there 
is another coupling mechanism that also creates orbital 
excitations. Namely, if the potential of the core-hole is 
not of Aig symmetry, it locally shakes up the occupied or- 
bitals, which can leave the Ti ions behind in an orbitally 
excited state. The corresponding scattering operator can 
be written in a general form as Oq = J^i Oi e^'i'^'. Due 
to the octahedral symmetry of the TiOg unit, the local 
scattering operator Oi can be decomposed into compo- 
nents {fK^)i transforming according to the rows k of the 
irreducible representations j of the octahedral group: 

= ^P,,.Af,^(/W). e^'i-^'. (2) 

The operators {fK^)i induce an orbital transition of cer- 
tain symmetry, the polarization factors Pj^k depend on 
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the experimental geometry. The factors Mj depend on 
the multiplet structure of the intermediate states and 
can be obtained from cluster calculations. We take only 
up to quadrupole components of {fK^)i into account (the 
higher multipoles contribute much less), and assume that 
all Mj are equal. The spectra of the operators (/i"'^)q are 
calculated within the superexchange model. Fig. 3b com- 
pares the result of the calculations with the experimental 
data. Here the energy scale riJsE was adjusted to 80 
meV. Both the lineshape and the scattering intensity of 
the observed inelastic peak are reproduced for all scatter- 
ing vectors. The only noticeable difference is the weak 
shoulder that appears in the calculations at about 120 
meV. This arises from single orbiton excitations, which 
in this case are not forbidden for the [100]-direction. 

In addition to the orbital superexchange model, we 
calculated the RIXS intensity in the opposite scenario of 
strong orbital-lattice coupling. The orbital excitations 
are now local crystal field transitions activated by the 
shakeup mechanism. We use the orbital ordering pat- 
tern and crystal field splittings (200 and 330 meV) of 
Ref. [13 to obtain the spectra shown in Fig. 3c. (Closely 
similar values were also obtained by other calculations 
[1, The crystal field transitions are localized and do 
not show any dispersion, but due to the g-dependence of 
the transition matrix element, the intensity ratio of the 
peaks can vary. With increasing q, a shift of the overall 
spectral weight to higher energies is thus seen, as in the 
experimental data. However, the calculated spectrum is 
comprised of two well-separated peaks and is hence in 
poor agreement with the RIXS data, which show a single 
broad feature. An ad-hoc reduction of the level splitting 
away from the predicted values [1,11, 1^ could improve 
the agreement of the lineshape, but this does not affect 
the g-dependence of the intensity, which exhibits a trend 
opposite to what is observed. 

In conclusion, we have determined the momentum de- 
pendence of the RIXS spectra of LaTiOs and YTiOa and 
compared the results to quantitative predictions of the 
superexchange and crystal- field models. We found that 
the former model yields better agreement with the ex- 
perimental data on YTiOa. On a qualitative level, we 
note that the g-dependent intensities of the orbiton ex- 
citations of YTiOa and LaTiOs (Fig. 2b) exhibit trends 
opposite to those of the corresponding spin excitations 
[2, y, as expected in superexchange models. Mixing of 
spin and orbital correlations by the superexchange inter- 
action also accounts naturally for the strong tempera- 
ture dependence of the RIXS intensity, which is difficult 
to explain in terms of localized crystal field excitations. 
An open problem is the spectral weight of single-orbiton 
excitations, which is predicted to be larger than exper- 
imentally observed. A possible explanation is that we 
have assumed the spins to be fully saturated, while in 
YTiOs the ordered moment is lower than l/is even at 



zero temperature, presumably as a consequence of joint 
spin-orbital fluctuations. Such fluctuations are also ex- 
pected to reduce the orbital order parameter, to which 
the single-orbiton weight is proportional. 
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